Individual and Multi Vortex Pinning in Systems with Periodic Pinning Arrays 
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We examine multi and individual vortex pinning in thin superconductors with periodic pinning 
arrays. For multi-vortex pinning we observe peaks in the critical current of equal magnitude at 
every matching field, while for individual vortex pinning we observe a sharp drop in the critical 
current after the first matching field in agreement with experiments. We examine the scaling of the 
critical current at commensurate and incommensurate fields for varied pinning strength and show 
that the depinning force at the incommensurate fields decreases faster than at the commensurate 
fields. 

PACS numbers: 74.60. Ge, 74.60. Jg 
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Vortex matter in superconductors with periodic pin- 
ning arrays has been studied for several years now start- 
ing with the pioneering work of Fiory et al. Q. Re- 
cently considerable renewed interest in this system has 
come about due to advancements in nano-Uthographic 
techniques in which arrays of holes |^||, defects or 
magnetic dots j^-p^ can be created where various pa- 
rameters such as the pinning strength, size, and period- 
icity can be carefully controlled. These experiments have 
observed interesting commensurability effects where the 
critical current shows a maximum and magnetoresistance 
a minimum when the number of vortices equals an inte- 
ger multiple of the number of pinning sites. In addition 
peaks in the critical current have also been observed at 
some fractional matching fields Many applications 
of superconductors require high critical currents and the 
nano-engineered periodic pinning systems may be able to 
provide optimal pinning. A key question that needs to 
be understood in order for optimal pinning arrays to be 
constructed is where and how the vortices are arranged 
at the varied matching fields. In particular, it is not clear 
whether there are multiple vortices located at individual 
pining sites above the first matching field or whether only 
one vortex sits at an individual pinning site, nor has it 
been shown how this ordering affects the observed com- 
mensuration effects. 

In magnetization experiments in anti-dot lattices it 
was observed that the critical current dropped abruptly 
beyond a specific matching field which depends on tem- 
perature. It was inferred that below this field the pin- 
ning sites were capturing multiple vortices that could be 



strongly pinned, and that beyond this field the additional 
vortices were located in the interstitial regions and were 
not pinned directly by the pinning sites. On the other 
hand imaging experiments using Lorentz microscopy 0| 
of vortices in periodic pinning arrays observed only one 
vortex being trapped per pinning site and showed that 
at the integer matching fields the vortices form highly or- 
dered crystals with the crystal symmetry depending on 
the particular matching field. Subsequent simulations 
of vortices in periodic pinning arrays where only one 
vortex can be trapped at a pinning site also produced 
the same types of vortex crystals observed in the Lorenz 
microscopy experiments at the matching fields [ p^ . In 
addition the simulations also showed that peaks in the 
critical current occur at those matching fields where the 
interstitial vortices can form a highly ordered structure; 
however, some matching fields did not show a peak in 
the critical current, and the magnitude of the commen- 
suration effect varied considerably at different matching 
fields. Recent experiments with small pinning sites where 
it is expected that only one vortex can be pinned at an 
individual pinning sites observed commensuration peaks 
in agreement with these simulations |4]. There is also ev- 
idence from decoration experiments [ 5 1 that multi- vortex 
pinning at individual pinning sites in periodic pinning ar- 
rays occurs. Other recent imaging experiments with Hall 
probe arrays have also observed multi-vortices at indi- 
vidual pinning sites and show that above a certain field, 
additional vortices sit in the interstitial regions Q . These 
experiments seem to indicate that when multi- vortex pin- 
ning occurs the vortices in the pinning site merge to form 
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a macro- vortex. 

Another open question in vortex pinning in periodic 
arrays is why the commensuration effects are only pro- 
nounced near T^- One proposal is that the bulk pinning 
increases at low T and washes out the effect of the pe- 
riodic pins; however, samples with periodic pinning al- 
ways show a much higher critical current than compar- 
ison samples without the periodic pinning even for low 
T, indicating that the periodic pinning is still effective at 
low T. 

In this work we present results from simulations to 
elucidate the response of the critical current with multi- 
vortex vs individual vortex pinning. We model the multi 
vortices as individual vortices with increased number of 
flux quanta. We find that with multi vortex pinning, 
peaks in the critical current of the same magnitude oc- 
cur at every matching field. For individual vortex pin- 
ning the critical current drops off after the first matching 
field. These results agree well with experimental results 
for vortices interacting with magnetic dots where with 
large dots multi- vortex pinning could be expected. We 
also examine the scaling of the depinning force at com- 
mensurate and incommensurate fields as a function of 
the pinning strength. For all pinning strengths examined 
we find that depinning force scales linearly with pinning 
strength at commensurate fields. For incommensurate 
fields the depinning force scales linearly for large pinning 
forces but shows a crossover to a faster than linear scaling 
with smaller pinning strengths. This result may explain 
why commensurability effects are most pronounced for 
high temperatures where the pinning is reduced. In ad- 
dition we also examine the scaling of the critical current 
at the second matching field for multi- vortex pinning and 
individual pinning and find that as the pinning strength 
is increased the critical current saturates when individual 
vortex pinning is present. 

We model vortices in a 2D superconductor interacting 
with a square pinning array. The overdamped equation 
of motion for a vortex i is 

f,;^^=f™-Hf7-M, = r,v, (1) 

Here, is the force from the other vortices, f^''^ is the 
force from the pinning, fc; is the applied driving force 
corresponding to a Lorentz force from an applied current, 
and 77 is the damping which we set to unity. We model the 
vortex-vortex interaction potential as Uv = ln(r) which is 
appropriate for thin film superconductors. We have also 
conducted simulations using the modified Bessel func- 
tion interaction potential if 0(7") which is appropriate for 
bulk superconductors. The total force on vortex i from 
the other vortices is )■ impose peri- 

odic boundary conditions in the x and y directions. For 
the long-range logarithmic interaction we use an exact 
and fast converging sum . The pinning is modeled as 
attractive parabolic wells of radius Vp. 



f7 = {fplrp){v^-ri)Q{\r,-vl\l\). (2) 

Here, 8 is the Heaviside step function, is the loca- 
tion of pinning site fc, /p is the maximum pinning force, 
and A is the penetration depth. The pinning is placed in 
a square array. The initial vortex positions are obtained 
from annealing from a high temperature where the vor- 
tices are in a molten state and gradually cooling to T = 0. 
After the vortex configurations are obtained the critical 
depinning force is determined by applying a slowly in- 
creasing uniform force. We simulate two models. In the 
first, the pinning sites are small and can capture only one 
vortex each, so for fields greater than the first matching 
field additional vortices are located between the pinning 
sites. In the second model, multiple vortices can be cap- 
tured by larger pinning sites, and we assume that the 
vortices in the pinning sites coalesce and form a single 
multiple-quantized vortex. The multi-vortices interact 
with other vortices as nln{r) where n is the quantization 
of the vortex. Beyond the first matching field we start 
the simulation by annealing with the number of multi- 
ple quantized vortices already fixed so that we do not 
actually model the merging of the vortices. The multi- 
quantized vortices will still feel the same pinning force 
as the individual vortices. This is a reasonable assump- 
tion if the core of the vortex is smaller than the pinning 
site, since the maximum pinning force is determined by 
the gradient of the potential energy of the pinning site, 
which should not be affected by the size of the vortex. 

In Fig. 1 we show the critical depinning force versus 
field for the case where multiple vortex pinning occurs 
(thick curve) and for the case of single vortex pinning 
(thin curve). For B/B^ < 1.0 the depinning force is the 
same for both pinning radii which is due to the fact that 
the maximum pinning force for the pinning sites are the 
same. A peak at B/B^ = 1/2 occurs which is due to the 
formation of an ordered state as seen in previous simu- 
lations and experiments. There is a peak in both curves 
at the the first matching field where the vortex lattice 
has the same symmetry as the pinning lattice. For fields 
greater than B/B^ — 1.0 the curves deviate, with the 
single vortex pinning curve dropping substantially while 
the multiple vortex pinning curve only decreases by a 
small amount. The height of the peak for the multi- 
vortex curve at B/ B^ — 2.0 and the higher matching 
fields is the same as that for B/B^ = 1.0 since in our 
model the multiple vortices experience the same pinning 
force as the singly quantized vortices. The multi- vortex 
curve also shows matching peaks at the fractional fields 
n/2. For the single vortex pinning curve the maximum 
critical depinning force at B / B^ = 2.0 is much lower than 
that for the multi-vortex case. The fractional matching 
peaks are also suppressed for the individual vortex pin- 
ning case. 

The vortex configurations where only one vortex is cap- 
tured per pinning site are the same as those observed in 
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previous experiments and simulations. The vortex con- 
figurations for the multi-vortex pinning case repeat for 
each matching field with a square overall vortex lattice 
at each matching field composed of multiply quantized 
vortices. 

In Fig. 2 we show the results of transport experiments 
for Nb films with square periodic arrays of magnetic Ni 
dots for two different sized dots. Here the period is 
a — 600nm and the upper curve has a dot diameter 
of d = 400 — nm and the lower curve d — 530 — nm 
and T = 0.98Tc. Additional information about the fab- 
rication of the arrays can be found in Ref[8]. The val- 
ues of p{H) for both sized dots is roughly the same for 
B/Bcf, = 1.0 with the p minima values at B /B^ = 1.0 be- 
ing the same. For B/B^ > 1.0 the d = 400 — nm pinning 
is greatly reduced as seen by the increase in p{H); how- 
ever, minima still occur at the higher matching fields. In 
contrast the p values for the d — 530 — nm show only a 
minimal increase with H and pronounced matching min- 
ima of approximately the same magnitude at each field 
are observed. This result agrees well with the simula- 
tions, suggesting that for the larger dots, multi-vortex 
pinning is occurring at each matching field, while for the 
smaller dots, only one vortex is being trapped at each 
dot and additional vortices are located in the interstitial 
regions. 

In Fig. 3 we show that commensurability effects be- 
come more pronounced for weaker pinning. We consider 
multi- vortex pinning for fp — 0.8 and 0.1. Fig. 3 shows 
fp/fp for the two pinning strengths. At the matching 
fields the depinning force equals the pinning force. At 
the incommensurate fields the relative depinning force is 
much lower for the weaker pinning sample. 

In order to quantify the behavior observed in Fig. 3, 
in Fig. 4(a) we examine the depinning force for the com- 
mensurate field B/Bcf, — 1.0 and the incommensurate 
field B/B^ = 0.64 for a series of simulations with var- 
ied fp. The depinning force for B/B^ = 1.0 decreases 
linearly with oc fp which is the individual pinning 
regime. Since the vortex lattice is perfectly symmetric at 
the matching field, the vortex-vortex interactions cancel 
so the depinning force is determined only by the pinning 
energy. For B/B^ = 0.64, where vortex- vortex inter- 
actions should be relevant, we observed that for large 
fp the depinning force again decreases linearly with fp] 
however, near fp = 0.2 there is a crossover to faster than 
linear decrease with fp oc /p'^^. We have found the same 
behavior for other incommensurate fields; thus, the dif- 
ference in the critical current between the commensurate 
and the incommensurate fields grows as fp is lowered. 
This behavior could account for the pronounced match- 
ing peaks only near Tc where the pinning is weak. We 
were not able to go to weaker pinning where one might 
expect a crossover to a collective pinning regime with 
f'p ^fl- 

To compare the effect of pinning strength on the criti- 



cal current for individual versus multiple vortex pinning 
we have conducted a series of simulations at B/B^ = 2.0 
for varied fp for the two cases. In Fig. 4(b) we show that 
fp for the multi-vortex pinning case increases linearly 
with fp. We have also found the same linear increase 
at the other matching fields for the multi- vortex pinning 
case. The depinning in the multi-vortex case occurs in 
one stage with the entire lattice depinning at once. For 
the individual vortex pinning case there is an initial in- 
crease in fp for low values of fp followed by a saturation. 
This saturation can be understood by considering that 
the pinning of the interstitial vortices is caused by the 
repulsion of the vortices at the pinning sites which is in- 
dependent of the pinning strength. The saturation point 
also marks the onset of a two stage depinning process 
where the interstitial vortices depin first followed by the 
vortices at the pinning sites at a higher drive. A two stage 
depinning transition at B/B^ = 2.0 has been observed in 
experiments for widely spaced hole arrays ||l^ . Below the 
saturation point the entire lattice depins simultaneously. 

In conclusion we have compared multiple and individ- 
ual vortex pinning in periodic pining arrays. We find that 
with multiple vortex pinning, peaks in the critical depin- 
ning force occur at every matching field with the same 
amplitude. For individual vortex pinning the depinning 
force drops markedly for B/B^ > 1.0; however, peaks 
at the matching field are still present. These results are 
in good agreement with transport measurements on pe- 
riodic magnetic dot arrays for different dot sizes with the 
smaller dot systems showing results consistent with indi- 
vidual pinning while the larger dot systems show results 
constant with multiple vortex pinning. We also show 
that the commensuration effects are more pronounced for 
weak pinning, with the critical depinning force decreasing 
linearly with decreasing pinning force while the depinning 
force at incommensurate fields decreases faster than lin- 
ear. For multi-vortex pinning the critical depinning force 
scales linearly with pinning force at all matching fields 
while for individual pinning there is a saturation effect of 
the critical depinning force. 
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FIG. 1. Critical depinning force as a function of vortex 
density for samples with square arrays of pinning sites. The 
upper curve (thick lines) is the depinning system with multi- 
ple vortex pinning for B/B^ > 1.0, Tp = 0.5 and lower curve 
is the depinning line for a system when only one vortex is 
captured at a pinning site. 
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FIG. 2. Experimental values of p{H) for Nb samples with 
square arrays of Ni dots for two different dot diameters 530 
nm (sohd hue) and 400 nm (dashed line) at T/Tc = 0.98. 
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FIG. 3. The normalized depinning force vs B/B^ for 
fp = 0.8 upper curve (thick lines) and /p = 0.1 lower curve 
(thin line) obtained from simulations. Multi-vortex pinning 
occurs for B / B^ > 1.0. Here the commensurability effects 
can be seen to be more pronounced for the weaker pinning 
sample. 
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FIG. 4. (a) The depinning force vs fp for B/B^ = 1.0 
(filled circles) and B/B^ = 0.64 (open circles). The depin- 
ning force for B/B^ goes linearly with /p as shown with the 
fit while the depinning force for B/B^ = 0.64 goes as /p "'^ 
as shown with the fit. For fp > 1.0 the depinning force for 
B/B^ — 0.64 crosses over to a linear behavior, (b) The de- 
pinning force for B/ B^ = 2.0 for multi- vortex pinning (filled 
circles) and for single vortex pinning (filled squares). For 
the single vortex pinning there is a saturation in the depin- 
ning force for fp > 0.3 while the multi-vortex pinning case 
increases linearly with fp. 
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